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Abstract 

The stellar (n,j) cross section of 40 Ca at kT=2h keV has been measured with a combination of 
the activation technique and accelerator mass spectrometry (AMS). This combination is required 
in cases where direct off-line counting of the produced activity is compromised by the long half-life 
and/or missing 7-ray transitions. The neutron activations were performed at the Karlsruhe Van 
de Graaff accelerator using the quasi-stellar neutron spectrum of kT=25 keV produced by the 
7 Li(p, n) 7 Be reaction. The subsequent AMS measurements have been carried out at the Vienna 
Environmental Research Accelerator (VERA) with a 3 MV tandem accelerator. The doubly-magic 
40 Ca is a bottle-neck isotope in incomplete silicon burning, and its neutron capture cross section 
determines the amount of leakage and impacting thus on the eventual production of iron group 
elements. Due to its high abundance 40 Ca can also play a secondary role as "neutron poison" for 
the s process. Previous determinations of this value at stellar energies were based on time-of- flight 
measurements. Our method uses an independent approach, and yields for the Maxwellian-averaged 
cross section at kT=30 keV a value of (cr)30fceV = 5-73±0.34 mbarn. 

PACS numbers: 25.40.Lw, 26.30.-k, 27.40,+z, 97.10.Cv 
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I. INTRODUCTION 



The doubly-magic 40 Ca is mainly produced in explosive burning processes during su- 
pernova explosions. In contrast to the "normal" hydrostatic burning phases the respective 
explosive burning is only ignited for a few seconds when the supernova shock front passes 
through the outer layers. The main fuels for the explosive burning stages are the so-called "a 
nuclei" consisting of 4 He clusters ( 12 C, 16 0, 20 Ne, and 28 Si). The ignition temperatures are 
between 1.9 GK (for carbon burning) up to 4 GK (for silicon burning). The time scales of 
neon and silicon burning depend only on the respective temperatures and not on the density 
since both are dominated by the photodisintegration reactions of 20 Ne(7,a) and 28 Si(7,a). 
Apart from its role in astrophysics, the (n, 7) cross section of 40 Ca represents an important 
test for the applicability of the Hauser-Feshbach statistical model in the limit of low level 
densities. 



A. Production of 40 Ca 

Woosley [lj] showed already in 1973 that a superposition of explosive oxygen burning, 
complete silicon burning with a-rich freeze-out, and incomplete silicon burning offers good 
fits to the solar abundances in the mass region 28 < A <62. 

In explosive oxygen burning a quasi-statistical equilibrium (QSE) in the mass range A= 

n n n 

28-45 is obtained at T>3.3 GK [2|, [3J, |4J). This temperature is too low to achieve full 
nuclear statistical equilibrium (NSE) but in zones where the temperature exceeds 4 GK 
minor production of iron group elements can occur. The main burning products are the a 
nuclei 28 Si, 32 S, 36 Ar, and 40 Ca. 

Explosive silicon burning occurs at T>4-5 GK and can be subdivided into incomplete 
burning, complete burning with normal freeze-out, and complete burning with a-rich freeze- 
out. Complete Si burning is attained for T>5 GK where full NSE is established and the 
iron group elements are produced. Complete Si burning with a-rich freeze-out occurs at 
low densities when the triple a-reaction is not fast enough to keep the helium abundance 
in equilibrium. Then traces of a nuclei remain which were not transformed into iron group 
elements. 

Incomplete silicon burning occurs at peak temperatures of T=4-5 GK when the tempera- 
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ture is not high enough for nuclear reactions to overcome the bottleneck at the proton-magic 
shell closure Z=20. The most abundant burning products are the same as for explosive oxy- 
gen burning, but partial leakage can produce iron group elements like 54 Fe and 56 Ni. Which 
of the three silicon burning scenarios actually takes place depends on the respective peak 
temperatures and densities during the passage of the supernova shock front. 

B. Neutron poisons in massive stars 

The neutrons for the s process produced in situ by the 13 C(a, n) 16 and 22 Ne(a, n) 25 Mg 
sources can be consumed by isotopes with high abundances Ai and/or large neutron capture 
cross sections. The most important "neutron poisons" are 16 and 25 Mg, the products of 
the neutron source reactions. These isotopes reduce the available neutrons per 56 Fe seed 
and thus the s-process efficiency Further neutron poisons are 17 0, 20 Ne, and 24 Mg. This 
list can also be extended to 40 Ca, since this doubly-magic isotope exhibits both, a rather 
high solar abundance which reflects the high abundance in the preceding generation of stars 
(see Table H]) , and one of the largest Maxwellian averaged cross sections among the 20 most 
abundant isotopes. 



TABLE I: Isotopic [f 


3] and solar abundances Ai (relative to Agi= 


10 6 ) [6] of the stable Ca isotopes. 


Isotope 


Isotopic abundance [%] 




Solar abundance Ai 


40 Ca 


96.941 (156) 




59200 


42 Ca 


0.647 (23) 




395 


43 Ca 


0.135 (10) 




82.5 


44 Ca 


2.086 (110) 




1275 


46 Ca 


0.004 (3) 




2.4 


48 Ca 


0.187 (21) 




114 


Sea 


100 




61068.9 



For the reasons outlined above, an accurate knowledge of the neutron capture cross section 
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of 40 Ca is of importance for the leakage through the N=Z=20 bottleneck and for the role 
as possible neutron poison. The previous recommended Maxwellian-averaged cross section 
(MACS) of (a) 30fceV = 6-7±0.7 mb at kT=30 keV |7( was taken directly from de L. Musgrove 
et al. [8( and is essentially based on the respective resonance parameters measured below 
-E' n =300 keV, with contributions from the total neutron cross section measurement of Singh 
et al. js], both using the time-of- flight (TOF) method. 

The remeasurement of the 40 Ca(n, 7) 41 Ca cross section with an independent method such 
as the activation technique plus following AMS could add confidence in these results. The 
very high sensitivity of the activation technique allows one to use very small samples so 
that neutron scattering corrections are completely negligible. A further advantage of the 
activation technique compared to the TOF method is that the direct capture component 
(DRC) is already included. In combination with accelerator mass spectrometry (AMS) this 
technique can be extended to hitherto inaccessible cases, where the product nuclei exhibit 
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Q). 



Direct atom 



a long half-life and no 7-ray transitions (e.g. 41 Ca, £1/2= 103000 a 
counting via AMS also avoids uncertainties associated with the 7-decay branching of the 
activation products. 

In Sec. [Til we describe the experimental techniques of the activation with a stellar neu- 
tron source and the subsequent counting of the produced 41 Ca with accelerator mass spec- 
trometry. Sec. II I II shows the data analysis, and the results are presented in Sec. IIVI The 
astrophysical conclusions of our results are given in Sec. B3 

II. EXPERIMENTAL TECHNIQUE 
A. Activation 

The activation was carried out at the Karlsruhe 3.7 MV Van de Graaff accelerator. Neu- 
trons were produced with the 7 Li(p, n) 7 Be source by bombarding 30 fim thick layers of 
metallic Li on a water-cooled Cu backing with protons of 1912 keV, 30 keV above the reac- 
tion threshold. The angle-integrated neutron spectrum imitates almost perfectly a Maxwell- 



Bo 



12 



tzmann distribution for kT = 25.0±0.5 keV with a maximum neutron energy of 106 keV 
. At this proton energy the neutrons are kinematically collimated in a forward cone 



with 120° opening angle. Neutron scattering through the Cu backing is negligible since the 
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transmission is about 98 % in the energy range of interest. To ensure homogeneous illumi- 
nation of the entire surface, the proton beam with a DC current of f«100 fj,A was wobbled 
across the Li target. Thus, the mean neutron intensity over the period of the activations 
was wl.O— 1.6xl0 9 s _1 at the position of the samples, which were placed in close geometry 
to the Li target. A 6 Li-glass monitor at 1 m distance from the neutron target was used 
to record the time- dependence of the neutron yield in intervals of 1 min as the Li target 
degrades during the irradiation. In this way the proper correction of the number of nuclei 
which decayed during the activation can be attained. This correction is negligible for very 
long half-lives, e.g. for 41 Ca, but becomes important for comparably short-lived isotopes 
and affects the determination of the 198 Au activity, since gold is used as the reference cross 
section for the neutron ffuence determination. 

A sample of calcium carbonate (152.6 mg CaCOs, corresponding to 8.926xl0 20 atoms 



40 Ca) with natural isotopic abundance (96.941% 40 Ca [13], see Tabled]) was used in this 
measurement. The sample material was pressed into a thin pellet 8 mm in diameter, enclosed 
in 15 (Ltm thick aluminium foil and sandwiched between two 30 /im thick gold foils of the 
same diameter. In this way the neutron ffuence in our experimental neutron distribution 



can be determined relative to the well-known capture cross section of 19 Au 12j. During 
the activation the gold foils and the Li targets were changed three times, subdividing the 
experiment into four single activations (see Table ILT1) . The net irradiation time was ~20 d, 
and the total neutron exposure was calculated using the activated gold foils to be 2.32 xlO 15 
neutrons. 

TABLE II: Overview of the activation parameters. "i a " is the irradiation time and "$t t" the total 
neutron exposure. For discussion of the uncertainties see Sec. IIII CI 



Activation 


t a [min] 


<S> tot [xl0 14 n] 




4075 


2.542 


ca-b 


10078 


7.085 


ca-c 


9909 


9.606 


ca-d 


4123 


4.008 


Total 


28185 


23.242 
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B. AMS measurement 



The number of 41 Ca nuclei produced in the neutron activation was measured offline 
applying the technique of accelerator mass spectrometry (AMS) at the VERA (Vienna 
Environmental Research Accelerator) facility. This facility is based on a 3-MV tandem 
accelerator equipped with a 40-sample caesium sputter source. Details on Ca measurements 



at VERA can be found in Refs. 
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15j . The sample material was sputtered in the ion source 



and the nuclides are analyzed according to their mass and energy. Calcium fluoride (CaF 2 ) 
was chosen as proper sputter material. It was shown to provide sufficiently large negative 
currents when extracting CaF 3 ~ ions, with the additional advantage that the production of 
the interfering isobar (KF 3 ~) is strongly reduced ltl [17I fis| . 

After the neutron activation, the entire Ca-containing pellet (CaCOs) was dissolved in 
nitric acid to ensure a uniform mixture of the irradiated sample material and, consequently, 
a cross-section value integrated over the full sample area is obtained. Several aliquots were 
used to prepare sputter targets by the following way: Under pH control, hydrofluoric acid 
(HF) was added to the solution and Ca precipitated as CaF 2 , separated via centrifugation, 
and then dried to get CaF 2 powder. After mixing with Cu powder this material was pressed 
into the sample holders and inserted into the ion source. Overall, 17 sputter cathodes and 
a similar amount of identical but non-irradiated blank samples were prepared for the AMS 
measurements. 

The Ca content was determined during eight measurement series applying two differ- 
ent particle detection setups reducing thus systematic uncertainties. AMS measures isotope 
ratios by determining count rate of the rare isotope ( 41 Ca) relative to the ion current of 
a reference isotope ( 40 Ca). The procedure applied at VERA is the following: Negatively 
charged ions are mass separated and injected into a tandem accelerator which runs at typi- 
cally 3 MV. At the terminal the particles have to pass a gas stripper canal where electrons are 
stripped off from the negative ions and parasitic molecular ions dissociated via collisions and 
Coulomb break-ups in the gas. The resulting positive ions are further accelerated. After the 
tandem accelerator the 4 + charge state is selected: i.e. particles having energies of 13.5 MeV 
and the right mass over charge ratio [Mj q) are further transported to the particle detection 
setup. The number of radionuclides is measured with an energy-sensitive detector, while 
the current of the stable ions is determined with Faraday cups. In case of Ca, fast switching 
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between stable isotopes and the radionuclide was applied at the injector magnet and both, 
40 Ca and 42 Ca currents, were registered. The isotope ratios of the unknown samples were 
compared with ratios of 41 Ca/ 40 Ca reference materials 19|. Blank samples, having negligible 
41 Ca concentrations, were used as controls for the quantification of the background level. 

Two different detection setups were used for these measurements, both with the goal 
to separate efficiently possible interferences in the signal of the radioisotope 41 Ca with the 
isobar 41 K. The first uses the so-called delta TOF-setup 20|, |2l). Here, the time-of-flight is 
measured for particles which pass a thick but very homogeneous, passive absorber foil made 
of silicon nitride (Silson Ltd., UK). K and Ca can be discriminated due to their different 
energy loss in the absorber, which is reflected in a different velocity when leaving the foil. The 
second method utilizes a recent development at ETH Zurich/Paul Scherrer Institut (PSI) 
of a compact type of ionization chamber which provides a high resolution in the energy 



measurement 



221 ]. Particles enter this detector through a sufficiently thick silicon nitride 



entrance foil, and isobars, dispersed in energy by the loss in the foil, are discriminated with 
a segmented anode via a AE — E measurement. 

For the measurement of the (n,7) cross section we expected an isotope ratio 41 Ca/Ca of 
~10 -11 , which is well above the machine background of 41 Ca/Ca at VERA and the level of 
blank and calcium material used as target for neutron activation. Both setups have been 
proven to provide a sufficient isobar suppression for our samples. 



III. DATA ANALYSIS AND RESULTS 



Determination of the neutron flux 



The measurement of the induced 198 Au activity after the irradiation was performed with a 
high purity germanium (HPGe) detector with a well defined measuring position at a distance 
of 76 mm surrounded by 10 cm lead shielding. The absolute efficiency for the 411.8 keV 
7-line of the decay into 198 Hg was determined with a set of reference sources and yielded 



e 7 =0.212%. 



The total amount of activated 198 Au nuclei, iVigg, at the end of the irradiation can be 
deduced from the number of events C(t m ) in the particular 7-ray line at 411.8 keV registered 



in the HPGe detector during the measuring time t m 23]. The factor t w corresponds to the 
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waiting time between the end of the irradiation and the start of the activity measurement: 

C(t T 

'7 J 7 ""Y ( J 

7 7 accounts for the relative 7 intensity per decay of the 411.8 keV transition (I 7 = 95.58±0.12 



^ 19S ~~ t — 1 — 7i m x t \ ~\~t ' (-*-) 

e~ i~ k-v (1 — e~ A tm ) e~ A tw 



% 



24j). For the measurement of the activated gold foils with the HPGe, the 7-ray self- 



absorption fc 7 has to be considered 23j]. This correction factor was 0.995 for all gold foils. 



The number of activated sample atoms N act is determined by 

N aet = N i {<T exp ){i)$totfb(i)- (2) 

In this equation, (a exp ) (i) is the cross section in the experimental neutron spectrum and <& tot 
the total neutron flux. The factor 

_ j l Q a (f>(t) e - A ^-*) dt 
fo 4>(t) dt 

accounts for the decay of activated nuclei during the irradiation time t a as well as for 
variations in the neutron flux. This factor is calculated from the neutron flux history recorded 
throughout the irradiation with the 6 Li glass detector at 1 m distance from the target. The 
reference value for the experimental 197 Au cross section in the quasi-stellar spectrum of the 
7 Li(p, n) 7 Be source is (<7 e:r p)( 197 Au)=586±8 mbarn [12 1. 

The time-integrated neutron flux Q tot = J <f>(t)dt seen by the sample (Table HIT) is deter- 
mined by averaging the neutron fluxes of the two gold foils enclosing the respective sample: 

tot N 197 (a exp )r 7 Au)i h {) 
B. Determination of the isotopic ratio via AMS 

Accelerator mass spectrometry determines the ratio of a radioactive isotope relative to 
one stable isotope of the same element. This measurement is carried out relative to reference 
standards with well-known isotopic ratios. In our case we have determined the ratio of 41 Ca 
versus 40 Ca, #=^- 

Fig. [1] shows the first and second energy loss signal from the ionization chamber measured 
with a blank sample and an irradiated sample with 41 Ca. The events from the quasi-stable 
isobar 41 K and the radioisotope 41 Ca can be clearly distinguished. With these settings, also 
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1250 1300 1350 1400 1450 

AE1 



FIG. 1: (Color online) Top/Middle: First (AE1) and second (AE2) energy loss signal from the 
anodes in the ionization chamber. The events from 41 Ca and its stable isobar 41 K can be dis- 
tinguished by comparing a blank and an irradiated sample with 41 Ca. In the top panel events 
from Se appear at higher energy losses which have the same M/q as 41 Ca. Bottom: Histo gram 
illustrating the separation between 41 K and 41 Ca. 

ions with the same M/q as 41 Ca can reach the detector, like 82 Se. However, as can be seen 
in the top part of Fig. [TJ these events are well separated from the A=A1 signals. 

The AMS measurements were repeated eight times. In the first six measurements the 
deltaTOF technique was applied, and in the remaining two the compact ionization chamber 
(see section II. B). The isotope ratio R=j^ from these measurements is (1.34±0.07)x 10~ n . 
The uncertainty is due to the 5.2% error (lex) from the AMS measurement. 

One can determine the experimental cross section (a exp ) for the 40 Ca(n, 7) reaction for our 
experimental neutron spectrum with the known total neutron exposure &tot= 2.324 xlO 15 
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from Table [TT] and Eq. [2j 



<^,X 4 °Ca) = ^ -J-. (5) 

Since ^ is the isotopic ratio R, and the correction factor f b for very long half-lives is equal 
to unity, this equation simplifies to 

(a exp )( 40 Ca) = R-^. (6) 

^tot 

From this equation we can deduce an experimental cross section of (a exp ){ 40 Ca) = 
5.77±0.34 mbarn (uncertainty including the total uncertainty from Table fTTTT) . 



C. Error analysis 

The experimental uncertainties are summarized in Table HTT1 Since stellar neutron capture 
cross section measurements are often carried out relative to gold [ij], the error of 1.4% [3] 
in the gold cross section cancels out in most astrophysical applications and was, therefore, 
not included in the present error analysis. An error of 2.0% was assumed to account for the 
uncertainty of the sample position (0.25 mm) relative to the gold foils during the activation 
that affects the neutron flux seen by the sample. The uncertainty of the HPGe efficiency 
calibration was derived from the accuracy of the calibration sources and of the calibration 
procedure. The uncertainty (1 cr) in the AMS measurements (5.2%) includes the statistical 
uncertainty, the reproducibility of the measurement and a systematic contribution from the 
measurement relative to reference materials (which is dominated by the uncertainty of the 
half-life of 41 Ca). The total uncertainty of 5.9% is the combination of the uncertainties from 
the neutron flux measurement (2.9%) and the AMS measurement (5.2%). 



IV. CROSS SECTIONS 



A. Evaluated cross section libraries 



The energy-dependent neutron cross sections from the evaluated cross section libraries 
JEFF-3.0A, JEFF-3.1, ENDF/B-VII.O, and JENDL-3.3 include experimental resonance pa- 



rameters adopted from the compilations of Mughabghab 
determined in a TOF measurement by de L. Musgrove [s 



25 



261 ] . These parameters were 



in the energy range between 2.5 
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TABLE III: Uncertainties for 197 Au and 40 Ca. 



Uncertainty (%) 



Source of uncertainty 


197 Au 


40 Ca 


Gold cross section 


1.4 a 


- 


Detector efficiency 


2.0 


- 


Neutron flux 


2.0 




Sample mass 


0.4 




7-ray intensity 


0.1 




Counting statistics 


0.2 




AMS uncertainty 




5.2 


Total uncertainty 


2.9 


5.9 6 



a Not included in the final uncertainty, see text. 

6 Including the 2.9% uncertainty from the flux measurement with the gold samples. 



and 300 keV with contributions from the total neutron cross section measurement of Singh 
et al. Q. The quoted uncertainty for most resonance parameters from de L. Musgrove fl 
is of the order of 10% and larger. 

The unresolved resonance region above 296 keV is based on different statistical model 
calculations and thus differs from evaluation to evaluation. For example, in JENDL-3.3 



this region originates from the code CASTHY 
average (cr) 30 fc e y 



271 ] and was normalized to the Maxwellian 
6.7 mbarn given in [tJ. The data sets of JEFF-3.1 and ENDF/B-VII.O 
are identical and use the recent nuclear model code TALYS [28], where the parameters have 
been adjusted to reproduce the existing experimental data. However, for the calculation of 
Maxwellian cross sections up to kT= 100 keV the contribution from the unresolved resonance 
region has only little influence and we used an average of all four data sets. 
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B. Maxwellian averaged cross sections 



In an astrophysical environment with temperature T interacting particles are quickly 
thermalized by collisions in the stellar plasma, and the neutron energy distribution can be 
described by a Maxwell-Boltzmann spectrum: 



$ = dN/dE n ~ ^JW n • e~ En/kT . (7) 

The experimental neutron spectrum of the 7 Li(p, n) 7 Be reaction simulates the energy depen- 
dence of the flux v •<& ~ E n -e~ En / kT with kT=25.0 ± 0.5 keV almost perfectly Q. However, 
the cutoff at E n = 106 keV and small deviations from the shape of the ideal Maxwellian spec- 
trum require a correction of the measured cross section (a exp ) for obtaining a true Maxwellian 
average, ((7)25 keV- This correction is determined by means of the energy- dependent cross 
sections from data libraries. 

Before calculating a Maxwellian averaged cross section, we determine a "normalization 
factor" NF which gives a direct comparison how good the agreement is between our mea- 
sured cross sections (a exp ) and the evaluated cross sections u(E) folded with the experimen- 



tal neutron spectrum of the 7 Li(p, n) 7 Be source 12j, (<j em /)=5.18 mbarn. Our experimen- 



tally determined cross section is (a exp )=5.77 mbarn, resulting in a normalization factor of 
NF=p2*L=l.ll. 

The respective cross section derived by folding the NON-SMOKER energy dependence 
0, [30I with our experimental neutron distribution would be (o"ns) = 12.48 mbarn, yielding 
a normalization factor of NF=0AQ. The Hauser-Feshbach model cannot be applied below 



a given minimum energy of kT '^40 keV (Tg ^0.48 GK) 30J due to the low level density of 
the double-magic nucleus 40 Ca which causes an overestimation of cross sections. Such low 
level densities occur also for very neutron- or proton-rich isotopes close to the driplines, low 
mass nuclei, and at low energies. 

Applying the normalization factor NF to the evaluated libraries affects also the corre- 
sponding thermal value and the strength of the individual resonances. However, since the 
uncertainty of the resonance parameters is, as mentioned above, 10% and larger, we decided 
to use this approach. As soon as more accurate resonance parameters become available for 



■10 



Ca, this data should supersede our way of extrapolating Maxwellian cross sections via 
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Eq.EJ 

{a)kT ~ A ' fE n e-^ dE n ■ (8) 

The factor 2 / ^/H comes from the normalization of the Maxwellian flux formula when using 
the most probable velocity instead of the mean thermal velocity. At kT=30 keV this yields 
(c)30fcev = 5.73±0.34 mbarn with our experimentally derived normalization factor NF— 1.11. 

We have also re-calculated the Maxwellian averaged cross sections using the resonance 
parameters given in ENDF/B-VII.O [31] with the R-matrix analysis code SAMMY 32], 
yielding (o") 30 fc e y= 5.13 mbarn at kT=30 keV. A calculation of the uncertainties was done 
independently by Monte Carlo sampling of the four most important resonance parameters 
which contribute with 85.3% to the Maxwellian average at kT=30 keV. Varying the res- 
onances at E n =10.83, 20.43, 42.12, and 52.57 keV within their given uncertainties yielded 
(c)30 fcev=5.04±0.30 mbarn. Varying the neutron widths T n of all resonances by 5% and the 
respective Gamma widths T 7 by 20% gives (17)30 fc e v = 5.16±0.31 mbarn, in perfect agreement 
with the SAMMY calculation. 

In Table IIVI these values are compared to the previously recommended values from Bao 
et al. which were adopted from Ref. 33] and originate directly from Ref. Qj. As can 
be seen in Table IIVI and also in the graphical comparison of the energy trends in Fig. [3], the 



differences between the values from [7, 
increase with increasing temperature. 



331 ] and the recalculated Maxwellian cross sections 



C. Comparison with previous values 



Table |V] compares in the upper part the extrapolated Maxwellian averaged cross sections 



at kT=30 keV from t 



grove Jsj and Singh [9|), Barrett [34|, and the Bao et al. compilation |7|, |33] (taken directly 



re evaluated library ENDF/B-VII.O 31] (using data 



rom de L. Mus- 



from [8]) with our activation cross section. In the lower part of Table IVl the resu 



four Hauser-Feshbach predictions are given. Incidentally, the value of Woosley 35] agrees 



ts from 



best with our experimental value of 5.73±0.34 mbarn, but also the recalculated value from 
the resonance parameters ((<r)30fcev = 5.16±0.31 mbarn) agrees within the uncertainties. The 
large discrepancies to recent NON-SMOKER 29, 30] and MOST (36, results indicate 
in the present case the limited reliability of the statistical model in view of the low level 
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TABLE IV: Maxwellian averaged cross sections (cr)kT (in mbarn) for 40 Ca calculated with the 
energy dependencies from evaluated libraries ("This work"). The second column shows the previous 

in the third column the results calculated 



recommended values based on the value IVont 

with resonance parameters in ENDF/B-VII.O 



3l| are given. 



kT 




(cr)fcT [mbarn] 




[keV] 


This work 


Ref. [7,8, 33J 


Ref. [31] 


5 


11.96 


11.3 


10.79 


10 


11.49 


12.1 


10.35 


15 


9.23 


10.3 


8.31 


20 


7.56 


8.6 


6.81 


25 


6.46 


7.5 


5.81 


30 


5.73 (34) 


6.7 (7) 


5.16 (31) 


40 


4.93 


5.8 


4.44 


50 


4.53 


5.4 


4.08 


60 


4.28 


5.2 


3.85 


80 


3.91 


4.9 


3.51 


100 


3.59 


4.9 


3.19 



density in the double- magic nucleus 40 Ca in this energy region. 
D. Statistical model and direct capture calculations 

Despite the fact that the Hauser-Feshbach model cannot (or only with large uncertain- 
ties) be applied below kT^AO keV we have performed statistical model calculations using 



the code "HFSM" with standard optical model parameters from Moldauer [381 . 1391 ] . The 
level density is based on the Fermi gas theory with pairing and shell corrections. The pa- 
rameters were fitted to reproduce the average level spacings of Mughabghab 25| . e.g. for 
s waves (D(/=0))=45±4 keV. The giant dipole resonance parameters for the calculation of 
the average radiative widths were £7=20.14 MeV with a width of 8.05 MeV and a peak 



15 



TABLE V: Comparison of Maxwellian averaged cross sections at kT=30 keV from previous exper- 
iments and theoretical predictions. Theoretical values are given without error bars. 





\ u /30keV 


Ratio to 




[mbarn] 


this work 


This work 


5 73 (34) 




ENDF/B-VII.O [31] 


5.16 (31) 


0.10 


Barrett [341 


10 (3) 


1 75+2-fI 

x - ' d -0.59 


Bao et al. [7, 33] 


6 7 (71 


■ L-J - ' —0.18 


Woosley [35] 


5.8 


1 ni+ - 06 


NON-SMOKER [29] 


13.0 


9 9 7 +0.14 
z - z ' -0.13 


MOST (2002) [36] 


24.1 


4 21+ - 27 

^• zi -0.24 


MOST (2005) [37] 


13.9 


2 43+ - 15 

z -^°-0.14 



cross section of 58.3 mbarn. These values yield an average radiative width of (r 7 ) s =193 
meV, which is a factor of ~8 lower than the 1.5±0.9 eV from de L. Musgrove {§]]. However, 
only with this value reasonable Maxwellian averaged cross sections can be obtained. This 
decrease is interesting since the large difference of almost one order of magnitude to the 
average p- and <i-wave widths ((r 7 ) p =360±90 meV and (r 7 )d=740±360 meV) discussed in 
Ref. 8j vanishes. 

Additional calculations were performed with the data set of ENDF/B-VII.O |3jj, based 
on the evaluation of Mughabghab 25|, [26] , and with the model code TALYS 28] for cross- 
checking. The results are consistent within a few percent for energies up to 10 keV, but 
differ by up to 20% for higher energies. TALYS calculations were repeated using a coupled 
channel code but this did not influence the neutron capture cross section. Thus the different 
energy dependencies from ENDF/B-VII.O and TALYS originated only from the different 
optical model potentials used. 

Direct radiative capture (DRC) calculations have been performed for s- and p-wave neu- 
tron transitions to several bound states in 41 Ca. Using the experimental thermal cross 
section of <Ti/ l =410±20 mbarn 26] and subtracting the contribution of 13.5 mbarn from the 
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tail of the s-wave resonances the estimated DRC cross section is 396.5 mbarn. The extrap- 
olation from thermal energies to the keV region leads to very small DRC contributions of 
0.4 mbarn and 5.4 /ibarn for s- and p-wave neutrons at kT=25 keV, respectively. 

The best estimate for the Maxwellian averaged cross sections in the keV region can thus be 
obtained directly from the resonance parameters given in 



25 



26j normalized to the present 



results. For neutron energies below E n =100 keV, Maxwellian cross sections are determined 
by single resonances, and the Hauser-Feshbach theory cannot be applied. The predicted 
Maxwellian averaged cross sections in Fig. [3] should therefore be taken with care and are 
only plotted for fcT>40 keV. 




tn 1CT 
in 
o 



10" 
1ff E 
10" 



-ENDF/B-VII.O 
- DRC (p wave) 

HFSM 

Talys 

■ NON-SMOKER x 0.462 




1 10 100 

neutron energy (keV) 



1000 



FIG. 2: (Color online) Energy-dependent cross sections <r(E) for 40 Ca from the evaluated data 
library ENDF/B-VII.O (black line) in comparison with the prediction from the Hauser-Feshbach 
models HFSM, NON-SMOKER (normalized by a factor of 0.462 to reproduce the experimental 
cross section), and TALYS. The red line in the lower part shows the DRC contribution from p wave 
neutrons. 



V. CONCLUSIONS 

We have measured the 40 Ca(n, 7) cross section via the activation technique in 
a stellar neutron distribution and derived a Maxwellian averaged cross section of 
(c")30fce^=5.73±0.34 mbarn at kT=30 keV. This value is, within the error bars, in good 
agreement with the previously recommended value of 6.7±0.7 mbarn [?], [3] taken directly 
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FIG. 3: (Color Online) Maxwellian averaged cross sections from this work compared with pre- 
viously recommended values 0, [3] and the values calculated from resonance parameters of 



ENDF/B-VII.O 3l[ (Table HV|) . The error band corresponds to the 5.9% uncertainty of our present 
measurement. Also plotted are the results obtained from the HF calculations HFSM and Talys, 
and the normalized NON-SMOKER predictions. 



from Ref. 8j, and higher than the value calculated from the resonance parameters 31] , 
(c")30fcev = 5.16±0.31 mbarn. With this result, 40 Ca still acts as a significant neutron poison 
in stellar environments, simply due to its high initial abundance. In the direct s-process 
flow 40 Ca is bypassed by its quasi-stable isobar 40 K. In thermally pulsing low mass AGB 



stars 



4(i ] Ca is therefore depleted by 13% in the s process. However, the further reaction 



flow via 41 Ca can reach equilibrium. Accordingly, (n, 7) reactions on 40 Ca contribute to 
the overall s-process efficiency in the Ca to Cr region. This holds also for the s process in 
massive stars, since the role as a neutron poison is not affected by (7, n) reactions, even at 
the high temperatures during C shell burning. 
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